Probability of the resistive state formation caused by absorption of a single-photon in 

current-carrying superconducting nano-strips 
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We have studied supercurrent-assisted formation of the resistive state in nano-structured Nb and 
NbN superconducting films after absorption of a single photon. In amorphous narrow NbN strips 
the probability of the resistive state formation has a pronounced spectral cut-off. The corresponding 
threshold photon energy decreases with the bias current. Analysis of the experimental data in the 
framework of the generalized hot-spot model suggests that the quantum yield for near-infrared 
photons increases faster than the photon energy. Relaxation of the resistive state depends on the 
photon energy making the phenomenon feasible for the development of energy resolving single- 
photon detectors. 
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I. INTRODUCTION 



A recently proposed hot-spot scenario of the resistive 
state formation in a thin, narrow, current-carrying super- 
conducting strip due to absorption of a single infrared 
photon* was subsequently suggested as a mechanism 
of single-photon detection^ in NbN microbridges. Al- 
though soon after realization such detectors found appli- 
cations, the first experimental evidence of the envisaged 
mechanism has been reported only a year ago*. The hot- 
spot modeli describes the formation and the following 
evolution of a normal hot-spot around the absorption site 
in an infinite homogeneous superconducting film where 
a single photon has been absorbed. The energy of the 
captured photon is quickly shared among a large num- 
ber of electrons via electron-electron and electron-phonon 
interaction leading to an avalanche-like multiplication 5 
of non-equilibrium quasiparticlcs (QP). The normal spot 
appears when the excess QP concentration becomes large 
enough in order to locally destroy superconductivity. At 
some stage the normal spot reaches the maximum size 
that is a trade-off between quasiparticle multiplication 
and their diffusion out of the absorption site. In a nar- 
row strip carrying a small sub-critical current, no resis- 
tive state is formed unless the maximum size of the nor- 
mal spot A n exceeds the strip width. If the spot-size 
is less than the strip-width, the resistive state may ap- 
pear due to a supplementary action of the bias current 
that should have a density j close to the critical cur- 
rent density jc- Once the normal spot becomes larger 
than the coherence length, the supercurrent is expelled 
from the spot. The effective cross-section available for 
the current flow decreases and the actual current density 
may exceed the critical current density. When this hap- 
pens, a normal domain spans the entire cross-section of 
the strip that causes a voltage pulse developing between 
the strip ends. The number of photons, which result 



in a voltage response, related to the number of photons 
crossing the geometrical area of the strip is a measure 
of the quantum efficiency (probability) of the resistive 
state formation. If A n is not sufficiently large in order 
to enhance the current density to its critical value, there 
is still a finite probability of the resistive state forma- 
tion due to multi-photon or fluctuation assisted events^. 
Since A n increases with the photon energy there should 
be a current dependent photon energy e that demarcates 
these two regimes of the resistive state formation. The 
threshold energy eo is defined by the simple geometric 
criteria A n (eo)/w = 1 — (j /jc) where w is the width of 
the strip. For photon energies e > eo but smaller than 
hfo (uq is the plasma frequency), the quantum efficiency 
should equal the probability of photon absorption and 
should not, therefore, depend on the photon energy. At 
e < eg the quantum efficiency is dominated by multi- 
photon and fluctuation assisted events, which generally 
have a smaller probability. 

So far the cut-off of the quantum efficiency at e « 
eo has not been experimentally observed, instead a 
monotonous decrease of the quantum efficiency with the 
photon wavelength has been reported^. Knee-like fea- 
tures have been found in the overall exponential depen- 
dence of the quantum efficiency on the bias current. Al- 
though they were attributed to the above discussed cut- 
off, the concluded size of the normal spot was far too 
large in comparison to model estimates 6 . The photon ab- 
sorption in a superconductor obeys Poisson statistics, in 
that the probability of absorbing simultaneously n pho- 
tons from a weak photon flux is proportional to the flux 
intensity to the power of n. Variation of the quantum ef- 
ficiency corresponding to n — 1,2 and 3 have indeed been 
founds in NbN microbridges. However, this observation 
would generally manifest quantum nature of any detec- 
tor and is not uniquely specific to any particular response 
mechanism. A time delay of the resistive state formation 



2 



in NbN meander lines has been recently reported^. The 
delay was attributed to the supplementary action of the 
bias current in the formation of the normal domain but 
the quantitative description implied the initial normal 
spot much larger than the hot-spot model predicts. 

In this work we report measurements of the quantum 
efficiency of the resistive state formation in narrow amor- 
phous strips made from Nb and NbN films. A reduced 
electron diffusivity and a reduced density of electronic 
states in NbN allowed us to observe the predicted cut-off 
of the quantum efficiency and to evaluate the intrinsic 
efficiency of the electron multiplication. We refine the 
hot-spot model and show that the resistive state forma- 
tion in the absence of fluctuations may be triggered by 
a single photon even though no normal spot appears. 
This approach eliminates contradictions within the ear- 
lier published experimental data and provides a quanti- 
tative explanation of how fluctuations enhance the quan- 
tum efficiency beyond the cut-off. We also present recent 
results showing that the resistive state created by a pho- 
ton "remembers" its energy, i.e. the effect has energy 
resolving capability. 



II. EXPERIMENT 
A. Sample preparation and measurement technique 

Experimental results were obtained with narrow me- 
ander lines made from either Nb or NbN thin films. Me- 
anders covering an area of 4 x 4 /im 2 were used in or- 
der to increase optical coupling with near infrared pho- 
tons. Contrary to the conventional technique, our modi- 
fied NbN films were prepared by dc magnetron sputtering 
of an Nb target in an Ar+N2 gas admixture that had a 
reduced partial pressure of N2. The deposition regime 
in combination with relatively high sputtering rates of 
about 1.2 nm/sec resulted in amorphous films with a ni- 
trogen content smaller than in the stoichiometric compo- 
sition. Films with a nominal thickness of 5 nm were de- 
posited on sapphire substrates kept at room temperature. 
Mostly because of the nitrogen deficiency, the supercon- 
ducting transition in 5 nm thick films occurred at w 6 
K that is about one third of the transition temperature 
of bulk NbN. Moreover, the transition temperature was 
partly suppressed due to the proximity effect between 
the central superconducting layer of the film and both 
the non-superconducting oxidized surface on one side and 
the film-substrate boundary on the other side. Films 
were patterned using electron-beam lithography and ion 
milling. Our Nb films had a thickness of 10 nm and 
were prepared on Si substrates by dc magnetron sputter- 
ing in pure Ar atmosphere. Films were patterned using 
electron-beam lithography and reactive ion etching. The 
line width was 84 nm and 98 nm and the normal state 
resistivity just above the superconducting transition was 
6 • 1CT 6 fl cm and 6.7 • 10~ 4 fl cm for Nb and NbN mean- 
ders, respectively. Due to a partial damage of strip edges 
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FIG. 1: SEM photo of a representative Nb meander on Si 
substrate. Visible fringes at the edges of the structure are 
formed by photoresist left for device protection. Inset shows 
a typical voltage pulse as recorded by the oscilloscope. 



by argon ions during the etching process, a portion of the 
meander line near the edges became normal. These nor- 
mal areas further decreased^ the transition temperature. 
For both materials the transition temperature Tq of the 
meander varied between 4 and 5 K. An image of a rep- 
resentative Nb meander made with a scanning electron 
microscope is shown in Fig. ^ 

In order to ensure a sufficiently high quality of our 
NbN structures we evaluated their superconducting pa- 
rameters relevant to the formation of a critical current 
state. Measuring the temperature variations of the sec- 
ond critical magnetic field B c2 near Tc, we concluded 
the diffusion coefficient D = 0.35 cm 2 sec -1 . This rel- 
atively low electron diffusivity assured that the mean- 
ders were in the diffusive, dirty limit. Extrapolating 
the linear temperature dependence of B c2 near Tc, we 
found B c2 (0) — 9.4 T and the zero-temperature coher- 
ence length £ = {^o/^B c2 (0)) 1/2 « 8.5 nm. The den- 
sity of electronic states Nq — 2.2 ■ 10 24 m~ 3 K _1 is then 
deduced from the Einstein's relation No = l/(e 2 p D) 
where p stands for the normal-state resistivity. Further 
assuming a classical energy gap Ao = 1.76 /csTc, we 
estimated the London penetration depth Ai(Q) = 1200 
nm. Since the density of states was an important pa- 
rameter for the evaluation of the size of the normal spot, 
we cross-checked it measuring the critical current den- 
sity as a function of temperature. Because the width 
of our structure was well below the London penetra- 
tion depth, the supercurrent density did not vary across 
the meander line. To derive the critical current den- 
sity from the measured values of the critical current we 
used the electrical cross-section of the meander. Due 
to the proximity effect at the lateral interfaces and the 
damaged edges the electrical cross-section was almost 
one halfS of the geometrical cross-section. The best fit 
of the measured temperature dependence with j(T) = 
meN^kBTcD/h) 1 ' 2 ^- (T/T C ) 2 )(1 - (T/Tc) 4 ) 1 / 2 
resulted in N = 2.7T0 24 m 3 K 1 in agreement with the 
density of states inferred from the film resistivity. Practi- 
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cal coincidence of the measured critical current with the 
depairing critical current evidences amorphous structure 
with a large concentration of defects that play a role of 
pinning centers. The granularity, if any, does not influ- 
ence the critical current and should not therefore effect 
diffusion of electrons. 

The substrate carrying the meander was thermally an- 
chored to the cold plate of a He 4 -bath cryostat and was 
operated at temperatures ranging from about 1.8 K to 
approximately 3.0 K. The meander was illuminated using 
an incandescent light source. To enable spectral measure- 
ments, the light was passed through a prism monochro- 
mator. The meander was voltage-biased through a volt- 
age divider mounted inside the cryostat. Estimated time 
constant of the bias was tb = 0.5 ns. The response to 
illumination was in the form of a random sequence of 
voltage pulses each of them signaling either absorption 
of a photon or a dark count event. Voltage pulses were 
amplified using broadband microwave amplifiers (noise 
temperature 6 K, band-pass from 0.1 GHz to 1.6 GHz) 
and then guided either to a 250 MHz bandwidth/1 GHz 
sampling rate oscilloscope or to a 200 MHz bandwidth 
voltage-level counter. Inset in Fig. ^ shows the resulting 
signal transient recorded with the oscilloscope. The mean 
photon count rate, that is proportional to the quantum 
efficiency, was normalized to the light intensity measured 
with either a silicon photodiode or a PbSe photodiode (at 
longer wavelengths). For dark count measurements the 
optical entrance to the cryostat was blocked completely. 

Relaxation of the resistive state, and consequently the 
voltage pulse across the meander, should last longer when 
initiated by a photon with a larger energy. Material pa- 
rameters of NbN, in particular electron-phonon interac- 
tion time, suggest that at 2 K the duration of the voltage 
pulse should be in the sub-nanosecond range. Bearing in 
mind the band-pass of our read-out electronic, we ex- 
pect it to integrate pulses in real time. The amplitude of 
the signal transient after the amplifier chain should then 
correlate with the energy of the absorbed photon. Vary- 
ing the discriminator level of the counter, we measured 
the statistical distribution of the signal amplitude and 
retrieved the mean value and the dispersion. 



B. Experimental data 

The normalized quantum efficiencies of the resistive 
state formation for Nb and NbN meanders are compared 
in Fig. [3 For both materials it was measured at 1.8 K us- 
ing a current / = 0.9/ci where Iq is the critical current 
at the operation temperature. The quantum efficiency 
for the Nb meander drops continuously with the increase 
of the wavelength showing no threshold features while the 
NbN meander demonstrates an almost constant efficiency 
of « 1% up to wavelengths that exceed 2 /im. Obviously, 
the threshold photon energy for Nb falls out of our ex- 
perimental range and should correspond to a wavelength 
less than 0.25 /xm. According to the model described 
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FIG. 2: Normalized quantum efficiency for Nb (closed sym- 
bols) and NbN (open symbols) meanders. 
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FIG. 3: Relative quantum efficiency for a representative NbN 
meander biased with currents 0.6 7c , 0.77 Ic , 0.8 7c and 
0.89 7c- (from bottom to top). Arrows mark wavelengths 
corresponding to the cut-off at each bias current. For conve- 
nience, the data sets are shifted arbitrarily along the vertical 
axis. 



in the next section, the threshold wavelength should be 
proportional to p(D) x l 2 (Aq cf) _1 . Given the difference 
in the resistivity, the film thickness and the electron dif- 
fusivity, more than an order of magnitude difference in 
the threshold photon energy between Nb and NbN me- 
anders is easily understood. In the sections that follow 
we concentrate exclusively on NbN meanders. 

A decrease of the current brings to light the cut-off 
of the quantum efficiency in the NbN meanders. This 
can be seen in Fig. |3| that shows the relative quantum 
efficiency of a representative NbN meander measured at 
2 K with different bias currents. For each current there is 
a plateau in the wavelength dependence of the quantum 
efficiency that is followed by a drop as the wavelength 
increases. The wavelength corresponding to the cut-off 
increases with the bias current. Alternatively, varying 
the current through a meander exposed to illumination 
at a fixed wavelength, one should also be able to observe 
the cut-off of the photon count rate. However, this way 
additional uncertainties are introduced in detecting the 
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FIG. 4: Normalized photon count rates for wavelengths 
0.8 fim, 1.4 /im and 2 /im (thick solid lines, from top to bot- 
tom) and dark count rate (dashed line) both measured for the 
NbN meander at « 0.5 Tc- The largest absolute count rate 
was 5 • 10 7 sec -1 . Arrows mark the currents corresponding to 
the cut-off at: a - 0.8 /im; b - 1.4 /im; c - 2 /im. Thin solid 
lines approximate the photon count rates beyond the cut-off. 

threshold energy. The current destroys the symmetry of 
the energy states of superconducting electron pairs. The 
pair states enter the energy gap and the minimal exci- 
tation energy for some electron pairs becomes smaller 
than the zero-current energy gap. The excitation en- 
ergy decreases with the bias current and becomes zero 
when the current density reaches the depairing critical 
current density. The decrease of the minimum excitation 
energy with the bias current smears the threshold fea- 
ture and introduces a current dependence of the photon 
count rate above the cut-off. Fig. 0] shows the current 
dependence of the photon count rate (thick solid lines) 
for different photon energies. Indeed, when the bias cur- 
rent is close to the critical value, all curves merge and 
the count rate slowly varies with the current. The curve 
corresponding to a particular wavelength starts to de- 
viate from this common dependence at a current that 
increases with the wavelength. We associate the devi- 
ation from the common dependence with the cut-off of 
the quantum efficiency. Threshold values of the current 
are marked by arrows. The current dependence of the 
photon count rate beyond the cut-off was approximated 
(thin solid curves) using the current dependence of the 
dark count rate shown as the dotted line in Fig. 21 The 
connection between the photon counts and dark counts 
will be discussed below. 



III. MODEL REFINEMENT 

Applying to our experimental data (Fig. |3J) the geo- 
metrical criteria of the hot-spot model, we find the maxi- 
mum normal-spot size of 8.5 nm created by a photon with 
the 1.2 /jm wavelength. An analytical approximation^ of 
the numerical solution of the two dimensional diffusion 
problem gives a close value of w 9 nm. According to 




FIG. 5: Schematics of the superconducting strip-line carrying 
a current I. Bold lines demarcate the smallest volume where 
a reduction of the superconducting electron pairs causes a 
change of their mean velocity. Grey cylinder depicts the cloud 
of nonequilibrium quasiparticles with the size A. 



this approximation, even smaller normal spot should ap- 
pear in conventional NbN films 2,3 because they have a 
larger density of states, a larger energy gap and a larger 
electron diffusivity. However, all those spots are compa- 
rable or less than the coherence length suggesting that 
no resistive state should appear in response to photons 
with larger wavelength. Contrary, our experimental data 
show that not only the resistive state appear but there 
is a well pronounced spectral cut-off corresponding to 
wavelengths larger than 1.2 /im. In order to eliminate 
this discrepancy, we will refine the model considering the 
effect of quasiparticles, which are not confined in the nor- 
mal spot. 

Assume a strip made of a superconducting film carries 
a supercurrent / < Ic as it is shown in Fig. [S] The film 
thickness d is chosen smaller than both the line width w 
and the coherence length £ assuring the applicability of a 
two-dimensional diffusion model for nonequilibrium elec- 
trons. The small film thickness results in a large magnetic 
penetration length that exceeds the strip width even at 
temperatures well below the superconducting transition 
temperature. In this geometry, the magnetic field of the 
supercurrent may be neglected and, in the absence of 
magnetic vortices, the local supercurrent density should 
be constant throughout the cross-section of the strip. 
However, since our strips are wider than the coherence 
length, one can not exclude the presence of vortices. If 
there are vortices they should obey a Kosterlitz-Thouless 
(KT) topological transition. Square resistance of our 
NbN films was always smaller than the resistance quan- 
tum h/e 2 resulting in the KT transition temperature 9 
only slightly below Tc- At typical operation tempera- 
tures « Tc/2 all vortices are then bundled into pairs 
and do not contribute to the resistance unless the bias 
current is very close to the critical current. Although 
bundled vortices modulate the supercurrent density, we 
will neglect this current distribution for a while and as- 
sociate the local supercurrent density j — e n s v s with 
the mean velocity v s and the mean density n s of paired 
electrons. The local current state may noticeably change 
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only if n s changes over a distance £ or larger along the 
current path, shorter perturbations are tunnelled by elec- 
tron pairs without energy dissipation. The smallest vol- 
ume of the strip relevant to a current change, the £-slab, 
is marked in Fig. [S] with bold lines. If n s in the £-slab 
decreases by an amount 8n s , the mean pair velocity in- 
creases and becomes 



(1) 



0.20 



- Sn, 



as required by charge flow conservation. The characteris- 
tic conversion time of the pair velocity practically equals 
the Ginzburg-Landau relaxation time rs h/A (A is the 
temperature dependent energy gap) and is small com- 
pared to the electron thermalization time r t h- Thus, the 
mean velocity in the slab instantaneously follows changes 
of the pair density. The £-slab switches into the normal 
state if the pair velocity exceeds the critical value v ac that 
corresponds to the critical current density jc = en s « sc 
in the absence of perturbations. 

The concentration of nonequilibrium electrons C (r, t) 
at a distance r from the photon absorption site evolves in 
time due to multiplication of electrons and their diffusion 



C(r,t) 



M(t) 
AirDdt 



exp{- 



ADt 



(2) 



where M(t) is the time dependent number of nonequi- 
librium electrons. At t — r t h nonequilibrium electrons 
have thermalized to the energy level A effectively be- 
coming quasiparticles. Their number reaches the max- 
imum value M(Tth) — (t/A where C < 1 is the effi- 
ciency of the QP multiplication. After the thermaliza- 
tion time Tth has elapsed and before subgap phonons ap- 
pear, the concentration of nonequilibrium quasiparticles 
locally equals the reduction of superconducting electrons 
Sn s . Assuming that well below the transition tempera- 
ture n s ~ N A, one can rewrite Q and find the small- 
est number of nonequilibrium QP that is sufficient for 
switching the £-slab into the normal state 



SN* = N A£wd(l -I/Ic)- 



(3) 



If the spread of the QP cloud is comparable to the coher- 
ence length, all quasiparticles contribute to the change of 
the condensate velocity and Sn s = M(r t h) /(w <i£). If the 
cloud is larger (see Fig. |5J) , only quasiparticles confined 
within the £-slab have to be taken into account. The 
size A{t) of the electron cloud is twice the radius of the 
spot that confines all nonequilibrium electrons. Equating 
M{t) to the integral of the electron concentration (Eq. 
@) over the cylinder with the radius A(t)/2 and the 
thickness d, we find the time dependent size of the elec- 
tron cloud A(t) = A[Dt ln(M(t))] 1 /2. if M{ Tth ) > 100 
and typical parameters of NbN films are considered, A(t) 
exceeds the coherence length already at an early thermal- 
ization stage t < T t h- For A 3> £ the absolute number SN 
of nonequilibrium electrons confined to the £-slab can be 
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FIG. 6: Multiplication efficiency of quasiparticles for differ- 
ent photon energies concluded from current (open symbols) 
and spectral (closed symbols) dependence of the photon count 
rate. 



evaluated analytically. Integrating the electron concen- 
tration (Eq. (J2J) over the slab results in 



SN = M(t) 



(4) 



The number of nonequilibrium electrons in the slab peaks 
at t « Tth when nonequilibrium electrons can already be 
treated as quasiparticles. The cut-off of the quantum 
efficiency occurs when the maximum number of nonequi- 
librium QP in the slab reaches SN * . The criteria becomes 
SN(Tth) = SN* that allows to determine the efficiency of 
the quasiparticle multiplication 
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Using material parameters of our films, we estimated an 
efficiency of w 0.1 in the near infrared range. 

The maximum size of the actual normal spot A n can be 
estimated solving the equation C(A n /2,T t h) = NqA and 
using the above value for the multiplication efficiency. In 
our films it turns out that A n < £ and according to the 
hot-spot model the formation of a resistive state should 
not be triggered. When applied to the conventional films 
with a stoichiometric composition^ this analysis even 
suggests that no normal spot at all is formed. Finally, 
it is worth noting the meaning of the size of the normal 
spot that is defined by the criteria A n /w = (1 — I j Ic)' 
it is the diameter of a cylinder that would confine all 
nonequilibrium quasiparticles at t — T t h if their concen- 
tration would have been uniform and equal AoA. 

Assuming that the values of the bias currents and 
wavelengths indicated by arrows in Figs. Eland El define 
the cut-off that satisfies the criteria SN = SN*, we cal- 
culated for each pair of current and wavelength the effi- 
ciency of QP multiplication. Fig.[S]shows the efficiency ( 
for different photon energies. The data points presented 
by open symbols were obtained from the current depen- 
dence of the photon count rate (Fig. BJ . They are less 
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accurate (see the discussion above) than the data ob- 
tained from the spectral dependence of the photon count 
rate. Therefore, we omitted the data point derived from 
the count rate of 0.8 /jm photons. The efficiency of the 
QP multiplication increases with the photon energy. Al- 
though the physical reason for the increase is not quite 
clear, we speculate that it might be due to a splitting 
of the electron cascade into the electron and phonon 
branches. The splitting occurs^ around the Debye en- 
ergy, which is slightly less than photon energies used in 
our experiment. 



IV. FLUCTUATIONS AND QUANTUM 
EFFICIENCY BEYOND THE CUT-OFF 

Besides the fundamental interest in mechanisms, which 
result in dark counts in two-dimensional films, fluctua- 
tions are important for sensor applications where small 
photon fluxes have to be detected. Depending on the 
cross-section dark counts in one-dimensional supercon- 
ductors may originate from the spontaneous formation 
of either classical^ or quantum^ phase-slips. When the 
cross-section confining the current path extends over the 
coherence length, dark counts can appear due to discrete 
fluctuations of the superconducting order paramete r 12 ! 13 
and/or number fluctuations in the gas of bundled mag- 
netic vortices^*. It has been observed in many experi- 
ments that the dark count rate strongly depends on the 
bias current. Since at temperatures well below Tc the 
density of the superconducting condensate does not vary 
much with the bias current^, we suggest that the current 
dependence of the dark count rate is connected only to 
the velocity of the superconducting electrons. In partic- 
ular, the velocity controls the minimal excitation energy 
of the electron pairs. Although variations of the pair 
density modify the current dependence of the dark count 
rate, experimental data showi^ that this modification is 
a week factor in comparison to the current dependence 
itself. According to Eq. JIJ, the bias current correspond- 
ing to the spectral cut-off at a fixed photon energy can be 
presented as 7* = 7c (1 — 5n s /n s ). For currents I < I* a 
photon counting event can be thought of as a dark count 
event associated with that portion of the superconduct- 
ing strip, where the density of electron pairs has been 
decreased due to photon absorption. The elevated pair 
velocity in this portion of the strip is determined by the 
relative decrease of the mean pair density Sn s /n s . In the 
absence of perturbations, pairs would have the same ve- 
locity at a larger current I\ = 1(1 — 5n s /n s ). Therefore, 
we suggest that the photon count rate at I < I* equals 
the dark count rate at the current 1% = I (Ic /!*)■ 

Thin solid lines in Fig. 21 simulate the photon count 
rates beyond the cut-off. They were obtained by scal- 
ing down the current dependence of the dark count rate 
(dotted line in Fig. HJ . The scaling factor was used as 
the only fitting parameter. Visually, the scaling factor 
appears as the normalized current corresponding to the 



cross-point of the fitting line and the level of the photon 
count rate at I — Ic (straight horizontal line in Fig. 0}. 
These values are marked with arrows and were used for 
the estimate of the QP multiplication efficiency, exclud- 
ing the one for a wavelength of 0.8 /im. The method holds 
as long as the number of electron pairs destroyed by the 
photon in the £-slab is not very large or, equivalently, the 
number of remaining pairs compares to their equilibrium 
number. If the pair density changes noticeably, the cur- 
rent dependence of the fluctuation probability will change 
like it does with temperature. We believe that break- 
ing of this condition did not allow us to satisfactorily 
fit the photon count rate for the most energetic quanta. 
An obvious spin-off of this approach is an understand- 
ing of the spectral dependence of the photon count rate 
beyond the cut-off. The spectral dependence should re- 
peat the temperature dependence of the dark count rate 
at a fixed bias current. Available experimental data 2 *^ 
qualitatively confirm this conclusion. 



V. ENERGY RESOLVING CAPABILITY 

The normal conducting domain created in the meander 
by the joint action of a photon and the bias current re- 
laxes via thermal phonons, which pass the excess energy 
to the substrate. Since the relaxation rate is fixed, the 
life-time of the resistive state in the meander and, cor- 
respondingly, the duration of the voltage pulse increases 
with the full excess energy deposited in the film by the 
perturbation. Thus the duration of the voltage pulse 
should correlate with the energy of the photon, which 
initiated the formation of the resistive state. In this sec- 
tion we show that the resistive state indeed memorizes 
the photon energy. 

A current -voltage curve of the NbN meander recorded 
in the voltage-bias regime at 2 K is shown in Fig. [7| 
After the critical current value is reached, the current 
drops signaling the formation of a normal domain. The 
point A corresponds to the smallest stable domain, which 
has a length of the order of the thermal healing length 
Lt = (Dte) 1 / 2 where te is the effective electron cooling 
time. When the voltage increases further, a plateau oc- 
curs at an almost constant current I m that maintains the 
domain in equilibrium. According to Ref^&, exactly at 
this current in a current biased meander the velocity of 
the domain walls would equal zero. The Stekly param- 
eter, which is a measure of the self-heating by the bias 
current, can be evaluated as a w (Ic/Im) 2 — 25. Us- 
ing the resistivity of our specimen and the dc resistance 
Rj 1.2 kil corresponding to point A, we estimate a ther- 
mal healing length Lt ~ 60 nm. With the diffusivity of 
3.5 • 10~ 5 m sec -1 this length implies an effective electron 
cooling time te « 100 ps. Using this electron cooling 
time, the measured critical current density 5.5T0 9 A m -2 
and the electron specific heat ce ~ 2 • 10 2 J m _3 K _1 , 
we independently estimated the value of the Stekly pa- 
rameter as it follows from the dynamic theorj^ of an 



7 



1M 




FIG. 7: Current-voltage characteristic of the NbN meander at 
2 K. The critical current Ic and the current 7 m equilibrating 
the walls of the normal domain are both marked by arrows. 
The area A (marked by the circle) corresponds to the smallest 
stable normal domain of a size approximately equal to the 
thermal healing length. 



electro-thermal domain a — pjci^E /ce(Tc — T). We 
found a ~ 34 in fair agreement with the value concluded 
from the I-V curve. The steady-state parameters de- 
liver a self-consistent description of the electro-thermal 
domain in our specimen ensuring reliability of the fol- 
lowing consideration. 

The Stekly parameter a 3> 1 implies that even when 
I <C Ic self-heating strongly influences the dynamics of 
the normal domain. The critical energy that initiates the 
thermal roll-off or quenching in a current-biased meander 
can be presented^ as 



Qo = L t wcIce(Tc -T)ai 2 ln(- 



;), (6) 



with 



T C (I)-T . I 
—— — — and i = — , 

T c -T V 



where Tq(I) denotes the temperature, for which the crit- 
ical current Ic(T) equalizes the bias current I. Assum- 
ing the standard temperature dependence of the critical 
current, we find for the NbN meander at I /Ic = 0.8 a 
critical energy Qo = 0.51 eV that is ~ 3 times smaller 
than the energy of a 0.8-/im photon. Depending on the 
time constant of the bias and on the strength of the neg- 
ative feedback, the quenching may partly occur also in 
the voltage-bias regime. To fully avoid quenching the 
time constant of the bias should be noticeably smaller 
than the electron cooling time. The rates of photon and 
dark counts at different biases are shown in Fig. 0] The 
highest count rate of » 5 • 10 7 sec -1 was measured at 
I ~ 0.98 Ic- Such count rate corresponds to a recovery 
time of w 20 ns that coincides with the full duration of 
the signal transient seen in Fig. U] We believe that the 
recovery time more than 10 times larger than the bias 




1 2 3 4 5 6 
Discriminator level (mV) 



FIG. 8: The rate of dark (closed symbols) and photon (open 
symbols) counts at different discriminator levels. Solid lines 
are best fits assuming normal statistical distributions of the 
pulse amplitudes. 




Photon energy 



FIG. 9: Spectrum of 0.8-fim photons recorded with the NbN 
meander at 2 K. 



time constant tb ~ 0.5 ns was most likely determined by 
the lower edge of the amplifier band-pass. 

At large bias currents very close to Ic the noise back- 
ground may be so high that dark counts outnumber pho- 
ton counts. Accordingly, uncertainties in the photon 
count rate increase. Therefore, we analyzed the statis- 
tical distribution of the voltage transients at a smaller 
bias I = 0.8/(7 when the rate of dark counts was at 
least an order of magnitude less than the rate of pho- 
ton counts. Fig. [8] shows the dark count rate and the 
photon count rate for 0.8-/zm photons measured at dif- 
ferent discriminator levels of the pulse counter. When 
the discriminator level is set low all events are counted. 
Increase of the discriminator level leaves more and more 
events uncounted until at a high enough level none is 
counted anymore. Assuming a normal distribution of 
the signal amplitude, we fit the integral of the distribu- 
tion to the experimental data. Best-fit curves are shown 
in Fig. [S] Two fitting parameters were used: the mean 
value and the dispersion of the normal distribution. We 
have found a best dispersion value 1.65 ± 0.03 mV for 
both dark and photon counts and mean values 1.7 mV 
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and 1.9 mV for dark and photon counts, respectively. In 
terms of the strength and statistical distribution, dark 
counts are identical to photon counts when photons have 
the critical energy Qo-The wavelength of a photon hav- 
ing this critical energy is 2.4 fim perfectly matching the 
longest wavelength we detected with this meander. Re- 
lating the mean values of the signal amplitudes to the 
energies of 0.8-/xm and 2.4-/mi photons we found a scal- 
ing factor of 0.14 mV/eV. Fig. |5] presents the spectral 
histogram for 0.8-/im photons retrieved with this scaling 
factor from the experimental data. The best-fit enve- 
lope shown by the solid line represent a normal distribu- 
tion with a dispersion of 6.5 eV. Although the effect has 
an energy resolving power, it is far too low for practi- 
cal purposes. Factors limiting the energy resolution are 
discussed below. 

We should first estimate the energy dissipated in the 
meander by the bias current. After a photon has released 
an energy larger than Q , the normal domain starts to 
grow. The propagation velocity of the domain walls, 
i.e. the superconductor-normal interfaces, determines the 
growth rate of the normal domain. The velocity decreases 
as the actual bias current drops following the domain 
growth. In our experiment the life-time of the domain 
roughly coincides with the time constant tb of the bias 
circuit. The size of the domain, and consequently the 
energy dissipated by the current, can be determined in- 
tegrating the Joule power with the particular current de- 
pendence of the interface velocity. Considering the sim- 
ple step-edge heat-dissipation model, it can be showni& 
that the interface velocity v increases linearly with the 
current: 

v = (a—)V 2 {i- im )e-y\ (7) 

TE 

where i m = I m /Ic- With this dependence the expression 
for the Joule energy dissipated within the time tb takes 
the form 



Using tb ~ 0.5 ns we found a current contribution of 
« 11 eV that is almost an order of magnitude larger 
than the energy of a 0.8-/^m wavelength photon. Given 
the linear dependence of the response duration on the 



total energy released in the meander, the ratio of the 
mean amplitudes of photon and dark counts (see Fig. |HJ 
correlates well with the energy contributed by the bias. 
Denoting by e\ and ti the energies of 0.8-/im and 2.4-/zm 
photon, respectively, we find (Eb + £i)/{Eb +£2) ~ 1-09 
in agreement with the amplitude ratio 1.1 concluded from 
the experiment. 

The dispersion of the life-time of the normal domain 
may originate from structural defects of the meander me- 
diated by the domain length. In NbN films defects are 
most likely boundaries between grains, which have a typ- 
ical size of ten nanometersii. Intuitively, the dispersion 
has a maximum value when the domain length is of the 
order of the grain size. Contribution of the film struc- 
ture to the measured dispersion seems plausible because 
the normal domain in our meanders covers only a few 
grains. The cooled amplifier also contributed to the mea- 
sured dispersion since its equivalent voltage noise was 
only slightly less than the amplitude of the voltage pulse 
delivered by the meander. With a SQUID amplifier that 
typically has a subkelvin noise temperature, the energy 
resolution should drastically improve. Finally, using a 
superconductor with a smaller energy gap one should be 
able to increase the conversion ratio of the photon en- 
ergy to the domain length. Given the same dispersion 
this will also lead to a better energy resolution. An esti- 
mate shows that the cumulative effect of all factors men- 
tioned should suffice for an energy resolution comparable 
to the resolution of conventional superconductor photon 
detectors. 

In summary, we observed the current dependent spec- 
tral cut-off in the appearance of the resistive state af- 
ter absorption of a single near infrared photon in nano- 
structured NbN superconducting films. We generalized 
the normal-spot model explicitly taking into account the 
diffusion profile of nonequilibrium quasiparticles. With 
this approach, we found a super-linear photon energy 
dependence of the quantum yield in NbN that was at- 
tributed to a splitting of the thermalization process into 
the phonon and the electron cascades around the De- 
bye energy. The resistive state was shown to memorize 
the energy of the photon, which had initiated it. Ulti- 
mate energy resolving capability of this effect should be 
comparable with the energy resolution of conventional 
detectors. 
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